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ABSTRACT 


The purpose of this thesis was to evaluate the accuracy of the existing Mini- 
Ranger network in Monterey Bay and to suggest means by which the network 
can be improved. This network consists of six Stations located around the Bay 
and installed by the Monterey Bay Aquarium Institute (MBARI) and the Naval 
Postgraduate School (NPS). 

In order to undertake this task, data were made available by MBARI which 
were collected on a cruise made by their vessel “Pt. Lobos”. Additional data were 
gathered on a second cruise by the vessel “Pt.Sur” of the NPS. 

The data analysis indicated network problems. An effort was made to iden- 
tify these problems and to compute various correctors. In addition, an equation 
has been derived which ena’vies use of Mini-Ranger data collected when signal 
Strengths are low. 

The estimation of the accuracies obtained from the network through the var- 
ious tests applied, and the conclusions drawn. can be used as a guide to future 


seis Of the system. 
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I. INTRODUCTION 


In 1988, the Monterey Bay Aquarium Research Institute (MBARI) and the 
Naval Postgraduate School (NPS) began work to create a permanent network of 
Mini-Ranger stations around Monterey Bay. The net was to provide precise po- 
siuoning for research vessels and other scientific operations. Initially MBARI 
was seeking coverage of as much of the bay as possible, especially the Monterev 
canyon feature, with absolute horizontal position acuracies on the order of a few 
meters. It was understood this would be difficult to achieve, given the equip- 
ment’s specified range (40 nm). This study was undertaken to determine what 
levels of accuracy and spacial coverage were being achieved, as the net was con- 
figured in early 1989, and also to suggest possible improvements. 

Positioning data from two cruises were examined, one from a “Point Lobos” 
@alise On 3) NiIAR S89, and a second from “Point Sur” cruise on 22 SEP 89 (the 
“Point Lobos” and “Point Sur” are vessels used by MBARI and NPS, respec- 
tively). Significant errors Were found in a high percentage of the position fixes 
from the “point Lobos” data set. [t was apparent that at least some of the prot- 
lem arose from the use of weak signals (low signal strength readings) and incom- 
pletely calibrated equipment. Consequently. additional measurements were made 
onshore over known length baselines in order to quantify the systematic and 
random error components of the range measurement under various signal 
Strength conditions.! These findings stimulate the direction of the rest of the 
Study: 

® What were the main sources of errors? 


e Is there any postprocessing that can correct the data and derive better posi- 
tions? 


® What is a realistic estimate of the accuracy achieved? 


Seexre tliere any chances for improvement? 


1 It should be noted that the “Point Lobos” Master station failed before the baseline testing 
was run. Thus it was not possible to confirm the error estimates denved from analysis of the cruise 
data with independent baseline measurements. 


Il. DESCRIPTION 


A. MINIF-FRANGER SYSTEM 

The Motorola Mini-Ranger system being used by MBARI and NPS is the 
Falcon 484. It is widely used for precise positioning in a Variety of nearshore 
marine applications like hydrographic surveying, dredging, oceanographic data 
collection. etc. Numerous excellent descriptions of the system are avaliable in the 
literature. 

The system operates on the principle of a pulse coded transpondei ime 
Master station (also refered as to as the Receiver- Transmitter) is located on the 
ship and ranges to a maximum of four Reference stations ashore. The maximum 
update rate is once per second. The Master and Reference stations transmit on 
different frequencies in the microwave band. 5570 and 5480 Mhz, respectively. 
Both range and signal strength data are output for each Reference station inter- 
rogated. 

Signal strengths are computed at the Master as a relative indication of the 
power level received from the Reference station. Values can varv from about 5 
to 99. The manufacturer warns that anv ranges with signal strength below 13 are 
“weak” signals, subject to random errors exceeding the 2 meters (1 o } Specifica- 
tion. Signal strength obviously falls off with increasing range. but can also drop 
due to anv degradation of propagation conditions between the ship and shore 
Station. Hence. weak signals can and do occur at any range. Monitoring of the 
received signal strength 1s critical to assessing the accuracy of all range data. 

The range is computed at the Master station based on the measured round- 
trip travel time for the pulse and a preset refractive index for microwave propa- 
gation in the atmosphere. Unless changed by the user, the receiver assumes a 
value of N =(n-1)x10°=320 as an indicator for the refractive index (n), which is 
the value used throughout this study. Output ranges also can be corrected for any 
constant (Systematic) error previously entered by the user. Because this constant 


error includes time delays at both ends of the lime; il 1s essential to havc a specie 


tl 


calibration of each Master/Reference station pair. The manufacturer notes that 
without calibration, systematic range errors of up to ten meters are possible. 

limissDOSSibie (0 liawemthe Teceiver correct ranges from “slope to horizontal” 
distances before output, based upon user provided elevations of the Master and 
Reference stations. During this study, both the MBARI and NPS equipment 
were configured to output “slope” ranges only, choosing to handle this correction 
as part of the position computation algorithm. 

The Reference stations are supplied with electric current of 22 to 32 Volts 
DC from a pair of common car batteries, or permanently from an AC supply 
through a converter. A possible insufficient current Supplv can cause errors in the 
displaved distance in the receiver [Mini-Ranger Operation Manual, p.JI-11]. 

The standard system has a maximum range of 37 km (20 n.m), with Refer- 
emee station antennas of 13 db. The maximum range can be extented to 73> km 
(40 n.m) by raising the antenna gain of the Reference stations to 19 db, if the 
Station’s elevation is sufficient (due to the dip of the horizon). The master station. 
which is on the vessel. has an omni directional antenna of 6 db. The discussion 
and results in this thesis are all given using extended (19 db) gain antennas on the 
shore stations. 

Another significant factor for a Reference station is the directivity of 1ts an- 
tenna. The given maximum range of the system corresponds to the main lobe di- 
rection and drops as one moves to the sides. Figure | shows the pattern of a 19 
db gain antenna of a shore station. We can see that the range drops under 30 
km at 40° off the center direction. This is significant for a station that uses its 


maximum range ability and thus makes the choise of antenna directivity critical. 


B. MONTEREY BAY NETWORK DESCRIPTION 

The Motorola Mini-Ranger network in Monterey Bav consists of 6 stations 
with high gain antennas (19 db). installed permanently on known positions 
around the Bay. The installation was done by both MBARI and NPS serving 
primarily the Biological and Oceanographic research of their vessels “Point 


Lobos” and “Point Sur” respectively. The positions at which the stations are in- 





Table 1. PRELIMINARY POSITIONS OF THE MINI-RANGER 
Sit iONS 


MAO a 
pee Uinicoordemnla Seals Eas 
Longitude(W) vation 
nates hactor 
deg min sec 
| 585260.161 
ier , 
aba ai. 4092490.284 | 227078 
rx 602947.672 
ats 99973 
i 9. 4085231.399 


36 49 44.407 609863.128 
eae ard 999749 
aa 121 46 04.912 4076611.345 Eo 


CLOT Gea TaN, 
4055915.264 


600434 .009 
IFT 24 
4051260.014 jr 


396669.490 
4051826.41 1 


999743 


99 FS 





Stalled are not 1deal due both to the permanent power supply problem and to the 
fact that some of the ideal station positions are on private properties. 

On board the “Point Lobos” the Mini-Ranger omni-directional antenna was 
meeied on the mast 9.1 meters above the water Ievel. Thus the maximum 
achievable range due to the line of sight effect was increased. 

The shores of the Bay and the nearshore areas are generally low sand dunes 
so some Stations don't have the proper altitude to reach their maximum range. 

Table 1 gives the preliminary positions [Schnebele. 1989] of the stations. the 
mie) erid coordinates and their scale factor. The horizontal datum is the 
American Datum of 1983 (NAD 83). The UTM projection used is for Zone 10 


(central meridian 123° W.) and the seale factor 0.9996. Any errors in the post- 


fA 


tions of these stations will propagate into the computed vessel position by an 
amount that depends on the geometry. If the errors are smaller than several 
decimeters then there is no significant effect. Revised positions are shown in Ta- 
ble 3, section IV-C. 

The stations are not ideally positioned and so may be moved to different po- 
sitions in the future to allow better coverage of the areas of interest. 


Station Doppler is not yet permanently installed for safety reasons. 


Hl, THEORETICAL APPROACH TO THE PROBLEM 


A. THE RANGING PROBLEM 

The technique used for computing the ship’s position using ranges from 
known reference stations will be derived and explained in the beginning so the 
reader will be able to follow and understand the data manipulation and results 
of this study. It’s well known that one way to handle the computations involving 
ranging data is to work on some appropriate projection (e.g the UTM projection). 
There are several techniques that can be used for the position derivation: 


e Averaging the coordinates resulting from all possible combinations of meas- 
ured ranges. 


e Position calculation using the two most accurate ranges (using of the rest as 
a check). 


® Graphical methods. 


In this thesis we use the variation of coordinates method which is described 
in Cross [1981]. For the position derivation, at least 2 ranges are necded (two 
unknown parameters in the plane solution X.Y). With more ranges one can eSti- 
mate the accuracy of the resulting position and also make statistical decisions 
about the quality of the data. 

Figure 2 illustrates the concept of ranging from the i-th station. Consider a 
Station S, in a coordinate svstem AY. The observed range, reduced to the plane, 
is OR, (not shown). Consider an assumed position (AP). The distance between 
the reference station and the assumed position is CR. The calculation of the 
values dX.dY, which represent the corrections that have to be applied on the AP, 
will give the desired true position (TP) after some iterations. 


As 1S obvious from the figure: 
(Aap = Aap — Asi 
OA as 1 ee, 


The key equation is the one for the calculated range: 


L 


dy 


a Motaes ae ae. Spay Be eee 
~~ 

ae 

Ts 


== — — wo = 





Figure 2. VARIATION OF COORDINATES METHOD 


2 oe wane 
(CR) = (A)ap + Oda 


In order to find the transformation from range coordinates to position coor- 


dinates it is necessary to find the derivative: 
2CRACR, = 2(A;) 4pdX 427 J par 
GCR — Waal. 


where 


Je 
Ve= i 
M, CR dx 
y (Vi )ap r 


The estimate of the true range then is equal to the observed range plus an 
simerevichn can be positive or negative, called the residual I’. It is also equal to 
the calculated range plus the increment dCR.. 


Miecwtite pance tien 1s OK,-+ 34 where , 1s the sesidual But its also 
CR, + dCkR,. Solving for dCR,, we take: 


ee Oi.) artes 
and finally: 
MdX + NdY = (OR — CR), + V; 


For the n ranges problem, in matrix form, the equation may be Stated as fol- 
lows: 


Ai=Bt+tT 


Where matrices A .B and & respectively are: 


AM, N, 
a We Me 
M, N, 
(OR — CR), 
R= (OR ~ CR), 
(OR 7 GK): 


The standard solution for this formulation using the least squares method 


(Bomford, 1980) gives: 
X,=(A'WAY'ATIVB 
V. = Ee —B 


Where the subscript a indicates that the value is approximate. The best sol- 
ution is obtained by iterating until X converges to less than some limit. 
W is the weight matrix. The standard daviation o, of a range measurement 


is an expresion of its accuracy. It is standard to weight the range measurements 





using —>- which ensures achieving the minimum variance solution [Hamilton. 
? 
1964]. 
Assuming that there is no correlation between measurements, the weight 


matrix becomes diagonal as follows: 


L/ot 0 , 0 

I/o3 0 

wal o Wee 2 | 
0 0 7 : Ifo; 


Motorola estimates the standard deviation (1 o ) of a single range imeasure- 
ment to be 2 meters. The individual user has to set his own estimated ( 6 ) for 
his equipment (see section IV-H). 

The unbiased estimate of the variance-covariance matrix of the adjusted po- 


sition 1s given by 


2, 
oat 2) 


This matrix includes the variances and the covariances in the coordinate sys- 


tem we use ( X,Y coordinates) and appeatseds 1olloxs. 


Sos oy a 1g 
= 2 (18) 
Es Oy 


Assuming that ¢, = 2m, then the combination of two ranges can result in po- 
sitional error of 7.8 m (drms) in the limiting case of 30° intersection angle. Drms 
is derived from the unbiased estimate of the variance covariance matrix as fol- 


lows: 


} 


Idris = Nic ae c*) 


This error of 7.8 m can be larger when the system is not calibrated. For un- 
calibrated systems the error of a single measurement can exceed 10 m [Motorola 
Mini-Ranger manual]. 

The quantity @¢ is the a posterior: variance of unit weight. It can be computed 


as follows: 


Where dof is the degree of freedom. or n-2 when solving for X,Y position co- 
ordinates using n ranges. Ideally. ¢; should equal unity. 

According to Cross [1981]. when the degrees of freedom are small it is dan- 
gerous to use a Single determination of ¢, to judge the quality of the position. It 
is better to use an average value of ¢, from the whole data set. 

But the big question is “how we can say that the value of ¢, is out of bounds? 
For this question a chi square one or two tailed test is the best indication. even 
Mioweh there are situations where it is not sufficient (Uotila, 1975]. If the test 
fails, this mav be an indication that one or more of the following error sources 
exist that cause the large value of a,. 


® Error in the mathematical model involved in the solution (Scale factors. 
NGinetice bmcex Used etc.) 


¢ Computational errors. 


e High correlations between observations or very poor geometry of the situ- 
ation. 


e Influence of the omitted higher order terms in Taylor series expansion which 
can be critical for higher order accuracy systems (primarily when the Master 
Station 1s far away from to the Reference staron): 


e Incorrect variance-covariance matrix of the observed quantities. Jt means 
that for a non well calibrated system, large values of o, are expected. 


e Blunders in observations arising from various factors in the measurement 
process (Multipath, unstable readings, etc). 

So if the statistical test used for ¢, shows that something is wrong an investi- 

gation is needed for the detection of the possible error source and its elimination 


irom tied ator 


B. PROPAGATION EFFECTS 

The line of sight limitation is one of the critical factors in the network instal- 
jation. The radius of curvature of the signal path 1s about 4 times greater than the 
Earth's curvature. so the heights of the Reference stations and the heisht oni 
Master station antenna govern the maximum range of the system, according to 


the formula, [Motorola Mini-Ranger Nlanual, 1981}: 


—_—. 


d=4.04(h, +h, ) 


Where /i, and /#,, are the heights mn the Reference and Master stations respectively 
and d the maximum range in kilometers. 
The maximum range of the system can be reduced due to the atmospigis 


conditions. Attenuation of the signal due to rain is given by [Casey. 1982]: 
Azar 


Where A is the attenuation of the signal in db, r is the rainfall rate in mm hr and 
a.b are frequency functions. For the Miniranger mean frequency (5.5 GHz), their 
numerical values are a= 1.48 x 10°? and b= 1.1469. 

Another formula is given for fog, based on the water concentration per cubic 


MVELCK: 


al 4.87x10 “Mfdb/km. 


Where M is the water content in gr/#i? and f is the frequency in Ghz. 

Table 2 gives the attenuation of the signal for rain and fog. at 5.5 Ghz which 
iS an average operating frequency for the Mini-ranger. In terms of distance, we 
have 8.0 km range reduction at 40 km for each db drop. So, for example. in a 
distance of 40 km with heavy rain conditions we have 1.2 db drop in the signal 
(11.6 km reduction in maximum range). This drop in signal power can introduce 
unexpected range errors when signal strength drops below the Critical Strength | 
wiineshold (CST). 


Table 2. ATTENUATION OF THE SIGNAL AT 3.5 MHZ 












Precipitation Attenuation Attenuation 
Rate (mm. hr) (db; km). (db/km) 


a 
a 
a 






Fog (gr/#2°) 






A possible problem that can appear in the survey area is the multipath phe- 


iemenon ol the sicnal, I]t cam cause constructive or destructive inteference at the 


recelver. There exist two types of multipath. The first, vertical multipath is 
caused by reflection of the signal on the sea surface. The destructive interference 
at the receiver due to vertical multipath 1s called “range hole”. The second one 
is the horizontal multipath due to reflection of the signal on walls or other sur- 
faces. The destructive interference of the horizontal multipath is seen as “unstable 
readings”. 

Range holes are verv important in networks that operate at line of sight and 
can be a major source of positional problems. The best way to avoid such a 
problem is to put an additional Mini-Ranger Master station at a different height 
Onmuilesvessel 

The accuracy of the range measurements 1s affected by the weather conditions 
(temperature, humidity, pressure). A fixed value of the refractive index (N) 1s in- 
put into the svstem prior to a survey and kept constant throughout the survey. 
The default refractive index (N =320) corresponds to atmospheric conditions of 
2026. 50% aia OS mp resmeChineln. 

Extreme atmospheric conditions. like high temperature along with either high 
relative humidity or dry atmosphere, can give rise to large scale errors in the 
measured distances. For this reason further investigation was done to estimate 


the magnitute of the atmosphericallv induced errors (See section IV-E). 


C. LOW SIGNAL STRENGTH 

several investigators [e.g Casev, 1982] have shown that ranse errorsssGae 
bias and standard deviation, tend to increase once the received signal strength 
drops below some critical level. This level, usually refered to as the Critical Signal 
Threshold (CST). is defined by the users accuracy requirements. Ivpicalcnoiees 
are between signal strengths of 8 to 18. The Motorola Falcon equipment itself 
flags signal strengths less than or equal to 13 as potentially suspect. 

The signal strength from the shore station 1s also a factor that can introduce 
large errors In the range measurements, but it is alwavs output by the svstem and 
can be used in subsequent data processing. When signal strength ts above the 
CST limit, the errors have a Gaussian distribution, (mean Value zero anemier 


tively small variance). But when the signal streneth is below the @S® jimmie 


o) 
_ 
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® 
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Figure 3. VARIATION OF MEAN ERROR WITH SIGNAL STRENGTHEI 


errors follow a logarithmic curve. That means the range error increases logarith- 
mically as one receives smaller values of signal strength [Casev. 1982]. 

Picune 3 derived by Casey, gives the CST at SS=S8 .. For particular applica- 
Wemescacl USer Can Set his own estimated value for CST in order to either reject 
the low signal strength meausurements, or to increase the standard deviation of 
the measurement. 

The SS of each signal depends strongly on the distance from the reference 
Station, (it drops with the inverse square of the distance). It can also vary due not 
only to multipath phenomenon, (as We explained earher), but also due to the an- 


gle to the main lobe direction and the atmospheric conditions. 


IV. TESTS ON THE MINI-RANGER NETWORK IN MONTEREY BAY 


A. MINI-RANGER CRUISE DATA 

On 31 March 1989 the Vessel Point Lobos, from MBARI, sailed for a trip in 
Monterey Bay to make trials and calibration on the Mini-Ranger stations. The 
Master Station was installed on the mast at height 10 m above the water. During 
this trip useful data were collected and position estimates calculated. This data 
was processed in a preliminary fashion in order to detect possible error sources. 

Figure 4 shows the routes of the vessel followed. A division of the Bay into 
smaller areas was used, (as shown on the figure), in order to subdivide the data 
for processing and error analvsis purposes. The division was based on the con- 
cept of nearshore and open b< Vv regions. 

Another trip was subsequently done for NPS research purposes on 22 Sep- 
tember 1989, and some more data collected as shown on the same figure. The 
master station for this cruise was different and the data w collected under differ- 
ent conditions of atmosphere, time and instrumentation. It was considered possi- 
ble that the positional results from these two data sets may vary in accuracy. 

First we examined the earlier data set of MBARI. which already had the 
positions and their statistics calculated. We were informed by MBARI that for 
the position derivation they used the method of variation of coordinates. No de- 
tails on the specific code in this software were avaliable. 

The first step was a general view on the residuals and the “standard dev1- 
ation” provided in MIBARI’s data of the whole trip. This showed that something 
did not work well when the vessel was nearshore , especially in the SW portion 
of the Bay. Figure 5 shows the areas of the greater and more variable standard 
deviation of the data set. In general as shown in Figure 5. areas E and Guia 
fewer problems. It was not obvious which station or stations combination was 
responsible for the big residuals. 

Another problem was the variation of signal strength and its unexpectedly 


high values. In manv portions of the data set, the signal strength was much higher 
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Ist DATA SET (MBARI) 
— — — — 2nd DATA SET (NPS) 





Figure 4. SHIP ROUTES.FOR THE TWO DATA SETS/AREA DIVISION 


SPORADIC BIG RESIDUALS 
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Figure 5. AREAS OF GREATER STANDARD DEVIATION 


than expected for such a range. The appearance of signal strengths greater than 
70 in distances more than 10000 meters is unusual. This is a strong indication of 
possible system malfunction. 

Because the data set did not show the expected accuracy, the following Series 
of tests and checks were performed: 


e A visit to the stations for inspection of directional limitations, antenna in- 
Stallation problems and other possible affecting factors. 


e A check of the station positions with geodetic survey techniques, (triangu- 
lation, traverse. or differential GPS), and a least squares adjustment of the 
Survey data with fixed geodetic contro] points around the Bav. 


e Check of the software used for the position derivation using the variation of 
coordinates algorithm with accurate grid scale and slope distance correction 
leemmidues, 


e A theoretical view of the geometrical accuracies with the use of various Sta- 
tion combinations in the Bay area. 


e A check on the effects of extreme meteorological conditions to check for the 
possibility of introducing range errors in excess of | meter. 


e A range hole prediction for the network. 


e An investigation into the possibility of correcting those ranges with signal 
Strengths below the CST. 


e A calibration test for determining biased or linearly varving errors and their 
Variability with signal strength above the CST. 
B. INSTALLATION CHECK AND ANTENNA DIRECTIVITIES 
On 20th of July. a visit to the stations TREVOR. WATS and PACKARD 
was done for visual] examination. check on the antenna directivity . possible sig- 
nal blocking from nearby obstructions and other possible affecting factors. The 


Stations HAYS and HANK were visited on the Ist of August and 25th of Sep- 





tember respectively. The alignments of the antennas were determined to approx- 
imately 2° using a small portable magnetic compass. All the magnetic bearings 
are given as true bearings. A short description is given below for each station and 
a map is presented (Figure 6) for clarification: 


e TREVOR: This station is installed on the roof of a house in Santa Cruz. The 
matenieais attached to the chinimey of the house and 1s oriented towards 
173°. The signal is blocked by a rocky hill from 215° to approximately 100? 


270° 


pop (7 mi) 


MAIN LOBE DIRECTIONS Monterey Bay Vicinity 
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Figure 6. OBSERVED ANTENNA DIRECTIVITIES OF MINI-FRANGER 
NETWORK IN MONTEREY BAY 
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(clockwise). There is no clear view of the Northern near shore Bay region but 
this area is of little interest. 


¢ PACKARD: Installed on the roof of a wood craft building. hooked on a thin 
mast. On private property between Watsonville and Moss Landing. Its an- 
tenna points to 250° and may be blocked by the Moss Landing electric plant 
Backs berwech 202 ame 205° Iihere were no other problems observed for this 
Station. 


e WATS: Installed on the roof of a commercial building near the beach at 
Watsonville. Its antenna points at 215° .The entire Bay area is unobstructed. 
The antenna 1s approximately 2 feet from the vertical wall of the building. 


See) So. Installed on the roof of Silas B. Hays hospital at Fort Ord with an 
Miimverupied Bay view and directivity of its antenna to 295° . The perfect 
view along with the height of the station makes it useful in some portions 
SW of Point Pinos . The signal is blocked trom the higher Pacific Grove hills 
mom 253° back to the South. 


¢ HANK: Installed on the roof of a house in Pacific Grove. While it has good 
melonit the Sifmaycannot reach the nearshore portion from Point Pinos to 
Monterev harbour due to high trees existing downhill. The directivity of its 
antenna is 345° and the signal mav been blocked by high trees from 005° to 
O35°. The station is offset from the position specified by the preliminary co- 
ordinates. This particular correction 1s given in the discussion in Section C 
below. 


See ore eCR: [his station was installed near NPS Beach Lab. but due to 
safety reasons was subsequenly removed. 

it was not possible to check the electric current supply of the stations, but a 

voltage stabilizer can help the system to be kept healthier and unstable outputs 


can be avoided in the future. 


C. POSITION CHECK 

Geographic coordinates for the station sites had been determined on the 
North American Datum of 1983 (NAD 83) by geodetic survey techniques, with 
ties to existing horizontal control points published by the National Geodetic Sur- 
vev [Schnebele, 1989]. Conventional ground survey and Global Positioning Sys- 
femieioPs) techniques were both used. [he survey data had been analyzed with 
a three-dimentional least-squares adjustment software package (GEOLAB, by 
me Osu. inc.). Absolute horizontal position accuracies estimated by the ad- 
justment were better than 8 cm (circular error. 95% confidence level). Additional 


work to confirm the offset positions of the Mini-Ranger units at PACKARD and 


ra 


WATS, has vet to be completed, but is not expected to shift their preliminary 
positions by more than a few centimeters at most. 

These station coordinates were converted to Universal Transverse Mercator 
grid coordinates (UTM Zone 10) for use in the vessel position computation algo- 
rithms. The UTM system has been adopted as the standard grid used by MBARI 
and NPS with this network. The specific Transverse Mereator Projection algo- 
rithm [Flovd, 1985] is accurate to within | cm for the entire zone. 

As described earlier, station HANK was offset from the preliminary position. 
The offset was measured by a simple magnetic compass for azimuth and a tape 
for distance and should have errors no greater than 5cm. The measured bearing 
was 066° true and the distance 41] inches. The two positions were at the same 
height, so the corrections in UTM are: dx =0.952, dv =0.424, dz=0.00 in meters. 

The results are Sumarized and the revised station positions provided in the 
Table 3. With the exeeption of station HANK, the revised positions of the 
Stations differ by Jess than 10 cm from the preliminary positions. It is concluded 


that the station positions are not the source of the observed inaccuracies. 


D. ESTIMATION AND COMPARISON OF DERIVED POSITIONS 
1. Development of position derivation algorithm 

The measured ranges were used to reeompute the ship’s positions by the 
method of variation of coordinates (Cross, 1981]. These reeomputed positions 
were then compared with those provided by MBARI. The comparisons showed 
that there is a discrepancy between the given positions and the reeomputed posi- 
tions on the first data set. which in some cases were significant. Because of these 
discrepancies, we ignored all the given MBARI positions and recomputed all their 
data. 

For this purpose. a program was written in Turbo Basie. This program 
uses the variation of coordinates method on the UTM projeetion suriace) Seam 
factors were eomputed for each line. The program is given in the Appendix | and 
can be followed by comments provided in the code. 

The measured distances from the svstem, for each case. are corrected for 


their slope with the assumption of a right triangle (Ggnoring the eurvature of the 


tl 
ty 


Table 3. REVISED POSITIONS OF MINI-RANGER 
STATIONS 


Latitude(N) 
Name Longitude(W) UTM coordinates | Elevation 
deg min sec 
Trevor SOS e a OO 585260.158 
| 122 02 31.536 4092490.287 
2. 
Wats ; 602947.636 
21 50 39. 4085231.429 


Packard 36 49 44.404 609863.147 
121 46 04.911 4076611.241 


36 38 33.82+ GUTOZIE29s 
121 47 45.894 4055915.205 


596670.450 
4051826.796 





Earth). This assumption inserts an error of 4 to 5 cm at the extreme ranges which 
is insignificant for Mini-ranger measurements in the areas covered by the net- 
work. but it makes the calculations simpler. For the complete correction one can 
look in chapter 15 of Laurila [1981]. 

The next correction that was made to the measured distances was 
for the scale factor of the UIM projection. The following formula computes the 


scale factor for the line between the Master and each Reference station: 


= 


m= imn,(1 +—) 
OR- 





tho 
ts 


Where m 1s the required scale factor, sm, is the scale at the central meridian 
(0.9996), R is the mean radius of the Earth (6375 km) and E, is the mean Easting 


value for each case. The mean Easting value is derived from the formula: 
D4 
lapel Tee eee 


Where £, and E; are the true Eastings (x; — 500000), of the Master and the Ref- 
erence stations respectively. The formula is a simplification of the one given by 
Bomford [1980, p193] and the error due to this is of millimeter order, which for 
our purposes 1s negligible. 

The outputs of the program are the least squares estimate for the position 
of the vessel, the variance covariance matrix L,. the residual vector (V) and the 
a posteriori variance of unit weight ¢?. Measurements were assumed independent 
so the diagonal weight matrix is treated as a vector in the calculations. 

The a posteriori variance of unit weight c? is the one that is used for as- 
sessing both the accuracy of the network and the ranging data. ¢? can be scaled 
as iS Well known, bv increasing the standard deviation of a single range meas- 
urement. The resulting ¢? is examined statistically with the chi squared, two 
tailed test, at a 95% confidence level. The null hypothesis. ¢? = |), was siesieg 
against the alternative 6741. When the null hypothesis is accepted. it is an in- 
dication. assuming no data blunders, that the system works within the specifica- 
tions provided by the manufacturers. When the null hypothesis is rejected several 
possible reasons for this exist. These were summarized in section III-A. 

It also should be noted. howeyer, that the presence of any unies@ieea 
systematic effects tends to invalidate this type of test. One potential source, not 
modelled in the position computation algorithm, involves the non-simultaneity of 
range measurements on a moving ship. In effect, the antenna is not stationary 
over a cycle of range measurements to all Reference stations. Given each range 
measurement takes about 55 milliseconds, four range measurements may be 
Spread over a period of about 0.15 seconds. At a nominal vessel speed of 6 knots. 


the possible error is about 0.45 cm. which is of little importance in this analysis. 


Additional error can be inserted due to mast motion. In rough sea condi- 
tions the error 1s big enough (more than 1 meter) and it has to be taken into ac- 
count by increasing the estimated standard deviation of the range measurements. 

2. First data set 

The results of the chi square test at a 95% confidence level using standard 
deviations of 2 and 3 meters on the measured ranges respectively, appear in Table 
4, 

The acceptable values of c? under these conditions are: 

For three ranges 0.00 < a2 < 5.024 

Pomnouienancesm0,029) = g- = 3.69 

The numbers given in the first two columns of Table 4 represent the per- 
centage of the number of positions for which the null hypothesis was rejected 
against the total number of the positions examined. The number provided does 
not include some abnormaly high a posteriori variances of unit weight which were 
@aleca Dy muUllipath effects, or low SS. These data were rejected as blunders. 
iiiem percentage appears in the last column. 

Eiecae be sechecxen en Gasol, themumber of rejections 1s clearly 
avem higher than the 5°o which one would normally expect. 

In the column ‘Difference’ we can see how many of the positions that 
were rejected when ¢ =2m. are accepted with 6 =3m. This 1s an indication ei- 
ther of poor calibration of the system or that the svstem does not perform io the 
Stameard Claimed by the manufacturer. I[t should be noted that the range data 
Maembecn Corrected for an estimate of SvStematic errors determined in a prior 
calibration. 

The results shown in Table 4 suggest the possibilitv that the master Sta- 
tion used mav have beeen poorly calibrated. For this reason, an examination of 
the distribution of the residuals for each reference station was undertaken. After 
examining the distribution of the range residuals for each reference station. it 
became clear that some stations were consistently biasing the position solutions. 
From the distribution of these residuals. the bias was estimated by comparing 


ie mean, median and most probable values. The mean residuals and the 
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Table 4. REJECTED PERCENTAGE WITHOUT REVISED 
CORRECTORS 


Rejected 
blunders 





adopted correctors are shown in Table 5. This simple averaging of observed res- 
iduals is not an 1deal procedure for estimating bias errors. [he least-squares 
process tends to redistribute an error in anv one range to errors (residuals) in all 
ranges. Recognizing this limitation. however, the procedure gives a useful esti- 
Mate Of WICOMmecteds bids crson 

The results for the corrected data are shown in Table 6 and can seuam 
rectly compared with Table 4. The improvement as we see 1s above 657m 

In the column ‘Difference’ we can see, in comparison with the one Or tor 
4. that the values are much lower, which js an indication that most of the call- 
bration uncertainty has been removed. 

Even after this procedure, some a posteriori Variances of unit We1eiime 
still high. especially in areas A and B. The possible reasons for this can be: 

e¢ Incomplete removal of the blunders. 


e Operation of some stations at low SS which in turn causes larger deviations 
from the normal values. 


e Incomplete removal of the bias errors described above. The technique used 
does not eliminate the entire error. but only estimates it. 


e For any one of a number of reasons the signal can be lost temporanie 
When it 1s regained 1t sometimes lacks stability for several readings. 


Table 5. MIEAN RESIDUAL AND VALUES USED FOR 
RANGE CORRECTORS 
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Table 6. REJECTED PERCENTAGE AND INIPROVEMENT WITH RE- 


VISED CORRECTORS 





e Higher value of the standard deviation of the system used than the one used 


for the weight matrix formation. 


3. Second data set 
The second data set was examined Statistically in the same way and the 
results appear in Table 7. We can See that the second data Set gives much better 
solutions. This was expected because it was known that the master station used 


for this data collection, was well calibrated. 


Table 7. REJECTED PERCENTAGE FOR THE SECOND 
DATA SET 


AREA Rejected percentage for o 





The previous data set. for reasons noted earlier, showed major problems 
when the ship was in areas A and B. The second data Set. however, Showeuegie 
sistent results in all areas and a low rejection percentage even when o =2m. 
These results suggest that the standard deviation of the range measurement 


should be between 2 and 3 meters. 


E. GEOMETRICAL ACCURACIES 

The position quality given from the Mini-ranger network in Monterey Bay 
depends strongly on the relative geometry of the ship and the reference stations. 
The geometric position quality can be represented in terms of error cllipses, con- 
fidence circles, or with several other appropriate techniques. When sevena 
Stations are used. computer techniques are necessary if detailed analvsis is to be 


undertaken. Such techniques are described in Perugini [1988]. 


The provided figures in Appendix 2 were taken from computer calculations 
using software associated with MBARI’s Navigation processor. The numbers 
provided represent the semi-major axis of the predicted error ellipse at 95% con- 
fidence level that is formed at the particular position, assuming the standard de- 
viation of a single range measurement is 3 meters. 

The best working combination appears in Figure 7 and is obtained using 
Stations TREVOR, PACKARD, HAYES and HANK. We can see that the ge- 
ometry of the position is very good with two combinations of only three stations, 
Pesewation combinations TREVOR, PACKARD, HANK and TREVOR, 
PACKARD, HAYES (see Appendix 2) 

AS Is easily seen, station TREVOR occupies a position that gives strong ge- 
ometry in the mid-Bay regions and its position 1s the most critical for the network. 
On the North side of Monterey Bay there are also other higher positions on which 
a permanent station could be installed with a much better view of the Bay. Mlost 
of them are in the University of California Santa Cruz campus and some 
arrangments have to be made. 

Im the Watsonville area the hills over Buena Vista Point have a good Bay 
view but the problem of a permanent electrical supply and Security has to be 
solved. Such a position would overcome the height problem of the WATS station 
while maintaining 1ts geometry. 

Stations PACKARD and HAYS are at the highest elevations tn their areas. 
There is no advantage to shifting their positions without deteriorating the geom- 
SIA 

ero DOPPLER docs not add much to the geometric quality of the net- 
work in the mid-Bav but can be used for positioning near Monterey harbor where 
signals from PACKARD and HANK are obstructed. It also could be shifted NW 
of Point Pinos on the lighthouse or on another high position to improve the ge- 
ometry in the Monterey canvon area. Further investigation on this station has 


to be done for the best choice of a permanent position. 


HANK 
DOP 


Monterey Bay Vicinity 





Figure 7. BEST WORKING COMBINATION OF THE MUINI-RANGER 
STATIONS 
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Station HANK gives very good geometry but its position has some limitations 
because of obstructions to its Bay view. Some station at the southern site of the 
Bay 1S necessary to maintain a good network geometry for the areas of interest. 

Generally we can say that the geometry of the network is not responsible for 
the big standard deviations observed in most regions around the Bay, but it be- 
comes more Important as one moves offshore near the Northen and Southern re- 


gions. 


F. ATMOSPHERIC CONDITION EFFECT OF THE AREA ON THE 
MICROWAVE PROPAGATION VELOCITY 

An Investigation was carried out to see the maximum possible error due to 
propagation velocity which could occur under extreme atmospheric Bay condi- 
tions as compared with mean Bay conditions. The atmospheric formulation is 
@esciived in Chapter 6 p.133 of Laurila [1981]. The effects of temperature (T), 
Pressure (P) and humidity (E) on the refractive index CN) are given trom the 
Essen formula: 
oe 10 


pa Qe 
ge 37.1922 





mere P and E are in mb and T ts the dry bulb temperature in A®?. But 


Fe (RANE, 
~~ 100 


Where RH is the relative humidity (%) and £, is the saturation vapor pressure, 
given in mb. at the dry bulb temperature. £, 1s a logarithmic function of temper- 
ouumerso the eficct on the refractive index depends on the temperature and the 
relative humidity as shown in Figure 8, for various relative humidities and con- 
Stant atmospheric pressure at 1013 mb. 

Daily weather logs. kept bv the Meteorology Department of the NPS, were 
searched for the maxima and minima of temperature and relative humidity in 
Monterey during 1988-89. This was done so as to compute the average variability 


of the refractive index and the possible error In distance that is introduced for the 
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Figure 8. VARIATION OF THE REFRACTIVE INDEX 


actual atmospheric conditions. The base for the comparison is the refractive index 
value (N = 320) that is used from the Mini-ranger for the range calculation. 

In Monterey Bay during the spring, temperature averages about 55 °F and 
humidity varies from 30° during the hotest days to 95% during the cold nights. 
During working hours, a relative humidity of 55% with temperatures 70 F° are 
usual. Extreme values of 95% RH and 55 °F are rare, but if they are accompa- 
nied with high atmospheric pressure, 1030 mb, then they give N = 343.6 which 
in turn corresponds to a distance error of 23.6 ppm. In a possible distance of 


4Q000 m the introduced error is about | meter. On the other hand, low atmo- 


a2 


spheric pressure of 1000 mb with 80 °F and 20% humidity, gives N = 288 and 
32 ppm error. which in distance terms means 1.28 m at 40000 m measured dis- 
tance. Average values of N for the spring during working hours are between 310 
and 330, which indicates that the assumed Refractive Index of 320 1s fair. 

During the summer the increased temperature is accompanied with a drop in 
the relative humidity so the result tends to be the same as in spring with an av- 
eeaee IN = 320 to 325. The extreme cases here also rarely introduce errors more 
than | meter at distances of 40000 m. 

The fall season is the one with a slightly different mean N and more vari- 
abilitv. Warm weather with high humidity is more hkely to occur. In the extreme 
cases, at maximum distance, the error could reach 2 meters. but it rarely occurs. 
iigeceaxcrase N = 330 is better to use. 

Winter shifts the refractive Index back to the lower values due to the lower 
temperatures of the atmosphere. which makes the extreme cases more rare. A 
Mean yalue of N= 315 is better for use. 

Generally the resulting crror should be less than a meter at extreme ranges 
and several decimeters on average, When extreme atmospheric conditions happen. 


Even though the error is small. the use of the seasonal refractive index is advised. 


G. RANGE HOLES 

Range holes occur due to the vertical multipath, which is caused from re- 
Meee ot the signal from the sea surface. This is common to all microwave svs- 
tems that use a line-of-sight propagation path. Figure 9 [Gilb-Weedon. 1976] 
illustrates the situation. 

Pre count the range holes from the first to the Nth. Here N represents the 
integer number of wavelengths difference between the direct and the reflected 
paths. The first range hole occurs at the longer distance from the shore station. 

Several factors affect the range hole characteristics: 


e The distance between the reference and the master stations. [he longer the 
distance. the wider the range hole zone. 


e The height of the stations. The higher the reference stations the fewer the 
range holes near them. 


REFERENCE 
STATION 


MEASURING DIRECT PATH 
STATION AEFLECTEO PATH * 


POINT 





Figure 9. RANGE HOLE OCCURANCE CONCEPT 


e The wavelength of the signal. This affects both the width of the range holes, 
and the distance from the reference station at which they occur. 


e The atmospheric conditions along the signal path cause loss of SS. Thus the 
reflected signal from the sea surface, due to additional losses from the re- 
flection. 1s too weak to interfere the direct one, when they both armme toma 
master station. 


e The sea surface roughness. This affects the tvpe of reflection which @eeins 
In a strongly wave modulated sea, the appearance of the range holes ts rare. 


e Changes in the sea level due to tides, winds etc., cause range holes to shift 
location horizontally. 

Small reflection angles cause more frequent range holes because the signal has 
lesser losses in comparison with the case when the reflection angle is bigger. The 
additional losses of SS, in the second case, weaken the reflected signal which can 
not then interfere with the direct one. 

lnpnealiie etic destructive interference is rarely total due to losses of the re- 
flected signal from the sea surface, but it can reduce the SS to below the seiisi- 
uvity of the receiver or below the CST in the case of Mini-ranger. This will result 


In signal loss, or worse, in a bigger standard deviation of the measurement. 


In the Monterey Bay network. the range hole predictions are shown in Table 
8 for each station. For the maximum range of each station we assume a ship an- 
tenna of 10 m. 

Range holes closer than the third occur rarely because of the bigger angle of 


incidence which attenuates more of the reflected signal. 


Table 8. RANGE HOLE PREDICTION FOR MONTEREY BAY 


Maxi- 
STATIONS mum 
range 


Range for the | Range for the | Range for the 
first hole second hole third hole 





Table 8 is based on the following formula [Gilb-Weedon. 1976], in which a 
Diame carth is assumed. KR represents the distance from the reference station to 


tieseamece hole center in meters. 





ly 


The distance of the maximum expected signal 1s: 
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R,, = 4h, Gk +1) 


Where /, and /, are the heights in meters of the referenee and the master stations 


Pe cet e mens tio wiavclenctl and k is an inteser (k = 1.2.3....): 
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Figure 10. RANGE HOLE SHIFT VS TIDAL HEIGHT 


In the Bay there ts a sea level variation due to tidal effects of 6 Toi: 
maximum. This effect tends to shift the range hole center towards the shore sta- 
tion at high water and away from the shore station at low water. The overall ef- 
fect of thts phenomenon for the stations appear tn Figure 10, for each range hole. 

Typical wind and wave conditions tn the Bay are such that ‘perfect’ signal 
reflections are not expected. Thus range hole problems should be rare and Itmited 
to conditions of calm sea accompanied by light winds. Range hole phenomena 


were not detected tn the positioning data studied. 
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H. CRITICAL SIGNAL THRESHOLD (CST) DERIVATION 

The relatively large distances and varving propagation conditions encount- 
ered in the Monterev Bay network result in many low signal strength readings. 
It 1s important that a method be found for using these readings if at all possible. 
Casev, 1982. has shown that low signal strength ranges can be used if corrected 
in postprocessing by an empirically derived relation between bias error, standard 
deviation, and signal strength. 

The recommended procedure to establish this relation is to perform repeated 
range measurements over a known length baseline, with a variable attenuator to 
ariilicially vary the recerved signal strength. These attenuators were not 
avaliable. but the effect was achieved bv ranging over relatively long baselines of 
Several different lengths. in order to obtain readings over a spread of low and 
high signal strengths. Twe sets of range readings were collected for different 
combinations of Master and reference stations. with signal strengths varving from 
about 5 to 25. 

For each Master:Reference station pair, ranges with signal strength over 15 
feweeayciaced tO estimate a mean range, free of low SS effects. Diffcrences be- 
tween this mean range and all low SS (< 15) ranges were computed, representing 
the additional error due to low SS. By this technique, error estimates from dif- 
ferent Master. Reference station pairs could be combined without concern about 
systematic differences between equipment pairs. The resulting range error esti- 
mates were grouped bv signal strength: 5-6, 7-8, 9-10 and so forth. Mean errors 
and standard deviations Were then computed for each signal strength group as 
shown in Table 9. As expected, both the mean error and standard deviation tend 
{0 increase as the signal strength decreases. 

The two data sets were compared to ensure repeatability of the results before 
deriving the empirical relation. In Table 9 we can See the results of the statistical 
tests for the two data sets. They are treated as samples from one population. but 
each with their own estimated means and standard deviations and examined with 


the Fisher - Behrens statistical test [Hamilton. 1964] at a 95°o confidence level. 
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as follows. The T test limit was taken from the Student’s t distribution table with 


degrees of freedom calculated by: 





3 2 
ae. 
nr i> ) 
OO 4 4 
0; 02 


— $$ -- oe 
ne (ny — 1) n3 (1 — |) 


The calculated value for the test is given by: 





Where ao? and of are the estimated variances, p, and pz, are the estimated mean 
errors, and fv, and #, are the number samples in the first and the second data sets 
respectively. 

We examine the hypothesis: 

H,iX, —X, = 0 The two data sets belong to the same population. 

Against the alternative 

Hy:X, — X, #0 The two data sets don't belong to the same population. Table 
9 shows the derived statistics from the two data sets and the limit of the T test in 
each ease. Everv calculated value, for the Fisher-Behrens statistical test. that 1s 
bigger than the T value means that the H, hypothesis is rejected so something else 
going on in the data sets. 

In the column ‘Calculated value’, the number 1s derived using the means and 
the standard deviation of the measu;ements. We can see that when SS ts 5-6 and 
17-18, the calculated value is bigger than the T test limit, which means that in 
these cases the null hypothesis IS Tejeerede 

The number given in the column ‘Calculated value with increased ¢ ‘ is a 
combination of the standard deviation of the long term variation in the instru- 


ment and its observed standard deviation. AS we can sce the result is a total 


38 


Table 9.) STATISTICAL TESTS FOR THE TWO SAMPLE DATA SETS 


Calcu- 


Ist data set Calcu- lated 

, 2nd data set 
(mean error / a) | . — lated value with 
(mean error/ @ ) 


value (u) | increased 





Statistical agreement of the data sets. We conclude, therefore, that they belong to 
the same population. 
The next step is the derivation of the equation that will permit the use of the 


lower SS more effectively. The equation derived was 


Pe eee 57206) 
lnx 
Where y represents the mean error in meters and x the low signal strength 
and the given numbers represent meters. The resulted y value represents the ob- 
served minus the calculated (C-O) range correction. 
Figure 13 shows this resulting best fitted curve. The correlation for the case 
is more than 0.999 which indicates that the resulting equation is verv close to 
perfectly modelling the data. This equation can be included in the software as a 


eemrection for the lower SS rather than rejecting them below the CST. 


oo 
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Figure Il. FITTED CURVE TO THE COMBINED DATA FOR TESTS 1 AND 
2 


Additionally, due to the increased standard deviation of the ranges at lower 
SS a higher value for standard deviation of a single measurement should be used 


In the weight matr.x formulation. 


I. VARIATION OF THE STANDARD DEVIATION WITH TIME 
Eigth sets of range measurements were collected over the known distance be- 
tween station HAYS and geodetic station RANGE 7. The choosen distance was 


short (2341 m) to eliminate the effect on distance of the variation of weatites 
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conditions. These 8 sets were collected on 8 different days and they have length 
of about 1 hour. with 3 second intervals between range measurements. Table 10 
gives the bias of these measurements with respect to the known distance and the 
standard deviation of a single measurement. The results from Table 10 show that 
the system has a standard deviation which is indicative of its precision rather 
than its accuracy. In addition, it has a variable bias. The total standard deviation 


for the system must reflect both these components. 


Table 10. STANDARD DEVIATION AND 
BIAS OF RANGE ERROR 





The biases shown in Table 10 were used Statistically to show that the stand- 
ard deviation of one measurement is 2.54 meters. which concurs with the value 
derived by evaluation of o2 estimates in section C of this chapter. The procedure 
followed for this derivation was a simple multiplication of the standard deviation 
of the mean bias found from the data set by the square root of the number of 
data sets. The user shouldbe aware that the standard deviation can be larger on 


board a vessel in a rough sea due to mast movements and the selectivity of the 
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Stations that the system makes(non simultaneous measurements to the reference 
Stations). The bias variation from day to dav makes the calibration of the 
Stations time consuming because a one day calibration will not remove all the 


bias: 


Table 11. MEAN DIFFERENCES AND STANDARD DEVIATIONS FOR 
STATIONS HAYS AND TREVOR 









Standard de- 
viation 
(Trevom 



















Mean differ- 
ence (Hays) 


Mean differ- 
ence (lI revor) 


Standard devi- 
ation (Hays) 





LO2 





J. VARIATION OF THE STANDARD DEVIATION WITH DISTANCE 

The second thing was the variation of the standard deviation with the dis- 
tance or with S.S, when it is above the CST. Several problems limited thtesaaim 
collection. like the unavailability of an attenuator. the poor network of known 
positions around the Bay and the low elevation of possible calibration sites. The 
data were sufficient only for the evaluation of station HAYS and partially of 
Station TREVOR. 

Table 11 shows the mean difference of the measured distances from (lemme 
distances and the standard deviations of each set of range measurements. We can 
see that there is variation of the mean error with distance in both stations but no 
indication exists that a function can be derived because the mean differences with 


the given standard deviations are not statistically different. Actually, a functional 


variation of the standard deviation with distance was not expected, when oper- 


ating above CST. 
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V. DISCUSSION AND CONCLUSIONS 


The Mini-ranger network in the Monterey Bay was installed for research 


purposes by MBARI and NPS and was examined from various aspects. Useful 


results were drawn that can help in the understanding of the operation of trans- 


ponder type systems in other areas, and that will make the present network more 


useful when operating at low SS. 


Figure 12 shows the overall coverage of the Mini-ranger with all the limita- 


tions taken into acount. except that weather conditions can attenuate the signal 


and in turn reduce the maximum range of the stations. For the maximum range 


derivation a height of the master station on the ship was assumed to be 10 meters. 


From this network and the data processed, the following conclusions can be 


drawn: 


For the calibration of such systems. if one wants the maximum accuracy that 
can be provided. a bias taken from single observations over known distances 
or with the baseline method may be not enough. It 18 better to use répeamed 
observations and a mean corrector derivation. 


The positions of the stations in the Bav were. with the exception of station 
HANK. determined to 4 sufficint level of accuracy such that they were not 
responsible for some of the large range standard deviations. 


The software of the user is important in the position derivation and can 
cause unexpected errors if not properly validated. 


The meteorological conditions in the Bay can cause attenuation of the signal 
and. in turn. reduce the maximum range ability. The range error iroineae 
delay of the signal due to meteorological conditions does not exceed tiem 
meters in the extreme conditions and 1s smaller than 0.2 meters undemaeme 
usual operational conditions. 


Range holes in the Bay can occur, but should be rare due to usual windy 
Bav conditions and the modulated surface from the swell. 


We have shown that ranges measured at low SS need not be rejected, but 
can be corrected using an appropriate error function. This function may vary 
from system to svstem. 


We show that the standard deviation can be considered to consist of two 
components. the unpredicted bias, which is smaller when the svstem 1s well 
calibrated and the measured standard deitation of the data. Generali 


4+ 





Figure 12. MINI-RANGER COVERAGE IN MONTEREY BAY 


o of one measurement under ideal conditions is a little greater than 2 meters. 
It is more proper to use o = 3 meters when the system operates at sea. 
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e There are no indications for variation of the standard deviation with dis- 
tance (or SS) above CTL. 
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APPENDIX A. COMPUTER PROGRAM LISTING 


REM PROGRAM Miknoce FILE: NICK3. BAS 
NICK KRIONERITIS Date: 8-8-89 
Modified 10-17-89 by KJS 
t ses'esfes'es estestesleslc Sesesleseveses' lewlewtou'e eta ese we sleves's le eseclevovocte le sle cle cle cteuledowectesleodoelesie sevevesevevedesodvedeves eveveves: 


REM This program calculates the optimum position of 
a ship using distances from known stations. It 
serves the solution for a miniranger network with 
accuracy pe er than 10 centimeter. - 

ve Yes'e'eve ses'e sevese ste slows verve sles! fouls ard s'e seslevevesle eslesievesleslosvevedeveveslesesloslevedoveseveves eslosie vertes'es'e slestesteste eve esis 
REM The subrutine SOLVERANGX was changed from CDR Kk. SCHNEBELE 
| program ''HYDROPLOT". 

Petede Veda Tete ese Tete eee ae ToT Te Te Te EET ETE TE TE TET TTS TS STE TS TS TS TSUN TE TS TS TS TS TS EOIN TS TCS TS TS TS TS TN TS TS ES ENR 
SouaChkh 9216 
SINCLUDE '"SOLRANGX. SUB" 
SINCLUDE "NICKMR. SUB" 
ON ERROR GOTO ERRORHANDLER 
Print “Enter the number of stations’ 
Lriput NG=" N% 
Lprint "Number of stations= " N% 


DECLARATION OF THE VARIABLES 

AP#(1): Approximate position of theship in meters (UTM) 
SP#(2): Station positions in 2D array in meters (UTM) 

R#( 1): Ranges observed on the receiver in meters 

LX#( 1): Stations Eastings for the scale factor affection 
"X#(1): False Eastings of the stations 

Mers¢ 1): Northings of the stations 

'  Z#(1): Elevation s of the stations from MSL 

BPy(1); Scale factor calculated for each case 

W!(1): Diagonal Weight matrix assuming indepedence 

among the measured ranges 

WODEZ(1): Index numbers of the 4 stations in each record 
| Paistay(2): UIM Coords (x,y,z) cf 6 stations possible to use 
1 
1 
1 
! 


Vectors for SOLVERANG2 but not used in this version 
CmOs!(1): Bias corrections for each of 6 stations 
Wets!(1): Weight factor for each of 6 stations 


patie P:( 1: 3) 

mimeo yYNAMIC SP#(1:N%,1:3) , ReCl: No - LG tN 

Pie wyYNANIC Saas Noam tl oe. 2 le 

memo YNANIC MC 1:N%) , SF#C1:N%) , BRes#(1:N%) pawl 12 NZ) 
DIM code%( 1: 4). Pals Carl 6 lees Ooms.) Gls 6s Woes tagle ©) 


PRINT "FREE ARRAY & STACK AT STARTUP", FRE(-1), FRE(-2) 
CALL GetStations(AllSta#(),CmOs! () ,Wets! ()) 

merint CmOs! (1); CmOs! (2); Cm0s! (3); Cm0s! (4); Cm0s 
O8LL SetFiles 

Print "Enter assumed position" 

Input "AX=",AP#(1) 
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Input "AY=" ,AP##(2) 
Input "AZ=" ,AP##(3) 
LPRINT ‘Assumed Position"; 
For i= iL <efeye Ss 
Print "Assumed position in UTM (X,Y,Z)="' AP##(i%) 
LPRINT USING "YHEHAHHEE. #"; APH#Ci%); 
Nex tea, 


Vo yededesodededededevedodcdedodedcsdcdededovedcdededededesdevedesevedededevedesedevedesededest 


REM Reena for Cae the subrutine SOLVERANG2* 


Bae lope be boi hae bac are kok he i eh ke eral Chai aie tir iariar ns eirlaefiel sctisc! <i we lagieekar irises rinc oki sls 
SFAC! =1. 0 ic for the subrutine for the 
SHGT#=AP#(3) corrections that have already done 
AP#(3)=0. do 
LM#=1. 1 _' Convergence Limit test (meters) 
Rite eee gee Mean ae of Earth in meters * 0.9996 
I dededevesidevededededdedededevededededededetedededekesededede sete 
REM: Declaration oF Riis Mat ea xe 
UOC sigan ge eels ar gg sere eres) vovs edeveve dese se ded ven 

For i=1 to N% 
W! Ci%)=. 25 


‘Diagonal matrix assuming 
independent measurements 


Next a. 

l eslededecevevevedscedededededesevescdevededesvedevedededesetedesededede 
Beep 2: CLS 
PRINT "Put Printer online in COMPRESSED PRINT Mode! "' 
PRINT ‘ Hit any key when ready to continue’ 
WHILE NOT INSTAT 

LOCATE 2.37 

PRIN 
WEND 


Print INKEYS 
LOOP TO READ AND PROCESS ALL DATA UNTIL END OF FILE HIT 


"WHILE NOT EOF(1) 
For kKZ=l) to 500 
Ale RG EU TEM gee, cane: oC) Ri HK) Oe) a eae 


1 ' -' ‘. You! te 
Carter beter ts seseses'es eve s'esle ses 


REM CO eas iene slope ane ¥ 
DEVE TERE TS TEE TS ETE Tee TET TE TET TEE ETS TE TE TE TE TE TEE TE TEE TS ETE TE 
For iv=1 to N% 
If (R#(i%) > 1.d0), Then 
DHGT#=Z7#( 1%) -SHGT# 
HGT#=Ri#( i% 2-DHGT# 2 
IF HGTi#< O. dO; Then 
PRINT ' Warning - at time ;secé&; 
PRINT " Hit any key to continue’ 
While Not INSTAT 
Delay 0.2 
Printevcn 
Wend 


tf ! ! 


DHGT exceeds Range’ 
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Print INKEYS 
Else 
R#(i%)=SQRC HGT} ) 
Ed TE 
Else 
R#(i%) =O. dO 
Edi lt 
Print Ranees corrected tor sllope—  R#(1%) 
Next i% 


a A Pa 


Vo gesesksecevesedevece vevedevevedeslesesvscdededs desesevecevesesesevesededesese 
Rerecalculation of the scale factors * 

Vo gedevedevevevedevededededesesevsdevevedsdededevedevededevervevesedevedededese 
AXP#=AP#(1)-500000 

For i=l to NZ 


wae =X#( i%)-500000 

My#( 1% )=AXP# 2 Ee eri) uae a) 2 
SFH( 4% I Om O96 O96 Mata.) / (Gokm: ) 
Print "Scale Pectors fe i UTM" SFC 1% 

iopint - Scale Pee eous for UTM="SF##(i% 


1 ers s': ale s'e od s'e se er Paula se a's s'es'es's oe we ee ae alee oe yark ae elen's aaale none Ke alea's se s'e Prey “) a's 


ll Correction for ein eae factors = 


1 tou! Ye aloaton? toatl te eloutaate 
fesleseslcsleveveses'e sesevevlevcvicslevesioveveveveseveses rity 


R#( 1% )=Ri#( ali «SEFC ee) 

Peimee Corrected rances= Rii(i% 
Lprint "Corrected ranges=" R##(i% 

Next i% 


1 ale ses'es'es'e s'es! en's sles‘es Sawlanteutsuion' sean sn'en'enten sn selon se eulon sn se se se salen sel selse se se seen oes 
4u du et a8 an CRG ARAL ARR AN ER EN AR AR ON RAL REN ER AN ER AN AS EN AR ARR IR ER EDAD ER AN AN C8 


. 


eslostes Ne slenies) sess sleeeslesle sje 


PoeerOormation of station position matrix * 
D devededevevedevedededevedesete tate sete dete cece de Tele NSAI AAI ASSN 
For i=l to N% 
#C iv, 1)=XH#C 1%) 
Pir i%,2)=Y¥i7( 1%) 
SP#Ci%,3)=0. dO 
Next i% 
1 
CALL SolveRangX(AP#() ,SP#() ,R#(),CmOs! (),N%,SFAC!, W!(), LM# ,BRes#(),_ 
Sx! CO S¥y' 46 Sty 1 Vanlance:) 
CALL beteool( seco ,AP#( ),code.() ,R#(),BRes#(),sxx! ,syy! ,sxy! ,variance! ) 


MeMeUSING — Adiifattty! Hake ge Hee es sec& APi#(1) AP#(2); 
For i% = 1 to 4 
LPRINT USING ™ sybyiet dtdeetet. dt tate Hs ~code%(i%) R#Ci%) BRes#(i%); 
mRext i% 
LPRINT USING " dys. dey. dt tHE. ot dE. «sxx! syy! sxy! variance! 
{ 
WEND 


Next k% 

CLOSE #1 

PROSE 2 

END 

ERRORHANDLER: 

PRINT "FREE AT ERROR" ,FRE(-1), FRE(-2) 

PRINT "AN ERROR TYPE", ERR, "HAS HAPPEN AT’, ERADR 
END 


LL Ga a) Sa atante a tenteu'sntoatsntaniaatsas'antsn'sntsutsn'se suse’ salsn wasn enlcn welsaisalontsatsutanes nan salansnisntsn'ensnten' sul sn' salon ons 
ESR EEA TOY na alla ale pe tar bale pr river kwrh ers ark ers ark rd rd ark rk er red vd ek ed ek ed ee kd ed rd ke ek 
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Y Sedesededcacdedededkicdedeadkvdedeacdedededededsdedcdkeaevedvacdede ved deseakiededeskdeskcsedkedeacdevecse ved 

INCLUDE "openfile. sub" 

tkdtkukikiek: CetStations Subroutine weerrevededdkdededs kjs: Oct 89 %* 
file: NICKMR. SUB 
SUB GetStations(Al1Sta#(2),Cm0s! (1) ,Wgts! (1)) 
LOCAL ies Se GCOde. 


CLS: PRINT "Insert Disk with File Si DAR ame) naveee. 
Print " Hit any key to continue when ready"; 
While Not INSTAT 
mecate 2. 3/ 
eGe uk quemmen ay 
Wend 
ee ita tele > 
GALL OPENFILEG@ 71, Bo Sh DAT eos 
For a, = 1 Cogs 
INPUT #1, code4,AllSta#(ik,1),AllStaf#(i4,2) ,AlISta#( 14,3) (CmOs!tat ee 
Wets! (i%),s$ 
Next i% 
CLOSE, 3-0 
END SUB 


MSebvetekitekarkt: SetFiles Subroutine  wewiwkkbikkietedsiedks ice Oct B89 wiswteevsctkes 
file: NICKMR. SUB 

t 

‘opens input data file of Miniranger values(as #1), and opens output file 
' for theveomputeduresniter tas a2) 


SUB SETFILES 
LOCAL s$,txt$,FiNames 


Clot =ss = S@Drrve: path name. ext) a 

PRINT "Enter MiniRanger data file for input” + s$; 
INPUT FiNameS 

CALL OPENFILE("1I'',1,FiName$ ,128) 


PRINT “Enter solutions file for output mrss, 

INPUT *txts 

CAL OPENED LGE@ 0 Anette en 

PRINT #2, "Solutions from file "; FiName$ 

END SUB 

l oedesecisstevedevedededs ReadNickMR Subroutine wWsitetdvedededessdeveevedededsss kjs: Oct B89 wisvedededs 
file: NICKMR. SUB 

specific to read mini-ranger files (from #1) from Krionertis thesis work 

with time in secs & 4 station numbers & 4 ranges. 


Pvt 
a ® 
1 
1 
1 
1 


SUB READNICKMR(sec&, code%( 1), rang#(1),X#(C1),Y#(C1),Z#(1) ,Al1Staj#(2)) 
LOCAL sSo a nuecode.() 
DIM DYNAMIC nucode%( 1: 4) 

{ 


DINE INPUD iaes> 

sec&=VAL(MIDS(s$ ,8,7)) 

For i, = ltor4 
nucode%(i%)=VALCMIDS(sS , 114+4*1%,4)) 
rang#(i%)=VAL(MIDS(s$ ,22+9*i% ,9)) 


"check nucode%(), if not same as code%() then get new station coordinates 
for the X#(),Y##(),Z#() arrays 


item code. 14) <= mucodes( 14) ).. bhen 
code4(i%) = nucode%(i%) 


X#(i%) = Al1lSta#( code%(i%) ,1) 
Y#(i%) = AllSta#(code%(i%) ,2) 
Pi, ) =e lista ( code.( 14) ,3) 
Ema 1 £ 
Bext i% 
END SUB 
 eveckvess AkekekkRiteSol Subrutine wWeewdededededesesdedededededededededed vedere kis: Oct BQw4vetere 
File: NICKMR. SUB 


1 
| Writes Miniranger Solution and error statistics to File #2 


SUB RiteSol(sec&, AP#(1),code%( 1) ,R#(1),BRes#(1),sxx! ,syy! ,sxy! ,vari! ) 
LOCAL ply 


PRINT #2, USING “dytsdédH+t AHPEBBHE. i HEHEHE. t's sec& AP#(1) AP#(2); 
em ee Fieeatatit. Ht ta. 
Fer i, = 1 to 4 

EAN =+'2, USING £5:. code%(i%) R#(i%) BRes#(ix); 


Next ie 

tt ; " ‘ 
PRINT #2, USING ~ #fyt. dHf afi. SHE dE. Eye HRYHE. dE 3 «<sxx! syy! sxy! vari! 
END SUB 
l Seveveveslevevevevesevedesevedesededes Soa Bol ae ar ae eee Sa ear ree ig tg ier eh rinaleise lari elem erin riage Grigg sist gear ie ar 


faa! Saw'a alo ulo nan -= 


t alee eae nee sesestes ’ en'e seve afsala velo e se "e rer esles 10 “' ele s'e s'e ss Sa wlealsalsntseiselseseselselse'selsese'se'sese'sealsa’se’ ae s'e a's estes wee salsa se ae 
hire ard wn ra veves @veven en riverier iw iwi wis rintiskterfarteriictiriertwts eveveveve ¢ 


‘Pile = = SOLRANGX. SUB 
tue ve ve ve ce te SOLVE RANGE POSITION (ver X) * * * SUBROUTINE * * * 


t 


‘Returns grid x,y in first 2 elements of apsni#() 
t 


"INPUTS: 
/apsnir(3) = approximate position as grid x,y,z 
spsns} = station positions iia Tatwavewas So1d x,y 42 


Pe@ista,,3) (array has nsta% rows by 3 columns for x,y,z) 


‘ORS#(nsta% o)= observed ranges to each station in spsns# row order 
rangcor! (nsta%) = C-O corrections for each range reading 


mista’ = number of stations/ranges 

‘scalefac! = grid scale factor (grid: true ratio) 

‘wets! (nsta%) = weights for each rangé (1/sigma 

'dxlim#} = convergence limit to stop iterating position 
solution (given in grid units - e.g. meters) 
SOUTPUTS; 

"apsn#(3) = x,y position returned in lst two elements 
(elevation of antenna, z, unchanged) 
"BRes#(nsta%) = residuals on observed ranges (ob-computed) 
"sxx! = variance in x (from Inv(ATWA)) 

'syy! = variance in y 

soa = covariance of x & y 

variance! = est of sigma-o (#tat FomMiRes-2))) 


SINCLUDE "LUDCMP. BAS" 
SINCLUDE 'LUBKSB. BAS" 
SUB SolveRangX(apsn#(1),spsns##(2) ,ORS#(1),rangcor! (1),nsta%,scalefac! , 


wets! (1),dxlim#,BRes#(1),sxx! ,syy! ,sxy! ,variance! ) 
LOCAL A#(),spsn#() ,ORa#, Ax#, Ay#, Bist ,ATWH#C) , ATWBH( ) ,ATWA#() , Indx%() , 
sum, iteration%,i%, j%,k%, YCol#() , ATWAInv#() ,xdx%(_), Btemp##() 
DIM spsn#( 1:3), ATWB#C1:2), ATWA#(1:2,1:2), Indx<( 12) xdx7, dee 
Btemp#(1: 4) 
i%=1 
For 3% = 1 to nsta% 
Tt ORS? (a7) > Se d0— chen 
xdx%(i%) = 4% "xref to non-zero ranges (range > 1m) 
UNGK a7 
long ol Ite 
Next 4% 
nranges% = i%-1 


DIM A#( 1: nranges%,1:2), ATW#(1: 2,1: nranges%) 
SolveAgain: 
For i% = 1 to nranges% 
ORaf#f=ORS#( xdx%(i%)) + rangcor! (xdx%(i%)) 
For j% = 1 to 3 
Spsil-( jn.) = Spsms7Cxaxcci.) ei 
Next 4% 
CALL RangeLOP( apsn#(),spsn#(),scalefac! , ORa#,Ax#,Ay#, Bi#) 
Avi(i%, 1)=Ax?# 
A#(i%,2)=Ay# 
BRes#( i%)=Bi# "Residual (o-c) on each range 
Next i. 


t 


' Compute ATW ¢2,nranges%|, & ATWB ¢2] 
1 


DOr joe =e leon 

sum=0. dQ 

For i% = 1 to nranges% 

ATWiH#( 3%, 1%) = A#(i%, 3%) “wets! (xdx%(i%)) 
sum## = ATW#( 4%,i%)*BRes#(i%)+sum# 
Next i% 

ATWB#( 4%) =sum# 

Next jie 


' Compute ATWA ¢2,2] 
t 


= 1 to nranges% 
sum} = ATW#H#(C4%,k%) *A#CK%, 1%) +sum#} 


ATWA#( 7%, 1%) =sum# 
Next 1%, }% 
t 


CALL LUDCMPCATWA#() ,2,Indx%( ) ,i%) 

CALL LUBKSBCATWA#(),2,Indx%( ) , ATWB#C ) ) 

sum!=0. d0 

BOY “14 =k) to 2 

apsni#(i%)=apsni?(i%)+ATWB#( i% ‘shifts position by dx,dy which has 
sum# = ATWB#( #% sums "been returned as ATWB() from LUEBKSB 


ay 
tv 


Next i% 
sum# = SQR(sum#) 
If sum# > dxlim#}, Then 
INCR iteration% 
If iteration% > 7, Then 
Beep 3 
Print "Position did not converge -- NO SOLUTION! "; 
Print Using "Delta Position = +4. ##i##meters”; sum# 
sxx! = 0.0: syy! = 0.0: sxy! = 0.0: var! = LOG10( sum) 


EAL SuUB 

Solel iise 

GoTo SolveAgain 
End ee 
' Compute variance-covariance values from inverse of ATWA 
sxx! = sigma-x Z 
syy! = sigma-y Z 
: sxy! = sigma-xy 
ERASE ATWi#, spsni "make some room in memory 


DIM YCoi#(1:2), ATWAInv#(1: 2,1: 2) 
EOujo— 1 to 2 


Orvis = 1 “to 2 Midembdey vector. Vol. to tind 
rio C17 = 02 dC ‘inverse using back-substitution 
Next i% 


YCo1#( j%)=1. do 
CALL LUBKSB( ATWA#(),2, Indx%() , YCol#()) 
Oise d5.—1 tO 2 
ATWAInv#(i%, j%)=YCol#(i%) 
Next 1% 
Next 4% 
CSNG( ATWAInv##(1,1)) 


" 

4 

=< 
| 


CSNG(C ATWAInv#(2,2)) 
SXy! CSNG( ATWAInv#(1,2)) 


Compute estimate of unit vafiigtm)/ msartek Res 


{ 
' Note that range residuals are in BRes#() 
! 


If nranges% = 2, Then 


variance! = 9999 ‘signifies no variance computable 
Else 
sum =0. d0 
For i% = 1 to nranges% 
sum# = Bres#(i%) 2*wets! (xdx%(it)) + sum 
Next i% 
variance! = CSNG(sum#/(nranges%-2)) 
End its 
If nranges®% < nsta%, Then "Set residuals into proper channel 


Peusi, = 1 to nstas 
Btemp?#(i%) = 0. d0 

Next ei. 

For j% = 1 to nranges% 
Btemp##(xdx%(j%)) = Bres#(j% 

Next 4% 

For iz = 1 to nsta‘o 


If nranges% = 2, Then Bres#(i%) = 0.d0 
Bres#(i%) = Btemp#(i%) 
Next i% 
End line 


END SUB 


Mae ke te ve ve od RANGELOP * * SUBROUTINE * *# * * *% ve WH Oe OH 
1 

‘Computes A matrix components Ax & Ay and grid range 

for input to multiple-LOP solution algorithm. 


‘Inputs: apsn#() = approx position x,y,z coordinates (grid) 
spsnj#() = station x,y,z coordinates 

scalefac!= grid scale factor (point or line) 

ORa/# = observed slant range 


' 
1 
f 


‘Output: Ax# & Ay# such that (Ax)dx+(Ay)dy = OR-CR 
OmCR# = Observed minus Computed range at apsn;}} 
"Note: Observed range (OR) has been corrected for elevation 
. differences and grid scale before comparing with CR. 
SUB RangeLOP( apsn##(1),spsn7#(1),scalefac! ,_ 

ORa# , Ax# , Ayi? , OmCR; ) 
LOCAL i%,CR#,dx##( ) 
DIM bee j1@ILeS)) 


FOr Mi. — sl tor 


dx#(i%) = apsn#(i%)-spsn#(i% 
Next.1% 
CR = SQRCdx#(1) 22Gx1;( Z)me2y) 


Ax## = dxi#(1)/CR# 
Ay# = dx#(2)/CR# 
dx##(3)=0. dO 
OmCR# = ORa#*scalefac! - CR# 
END SUB 


ss oa ale ale ntaatees sf as toate ato efactentacs os a's a welean as Pry ao ale A 
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' Sri rT a OS Ee 
Vesvsdetedetdotdevelioioiniioonsonnvoonvonionininiccoinhoncobinnnhonbbdolselvevsvedrlslsrdes: 
t a 
file: OPENFILE. SUB 

1 
q ¢ « « « o 

Opens a disk file using TBasic commands. Preforms error checking and 
t : 

allows recovery from common mistakes. 
! 


‘INPUTS: modes$ = One character designation for filemode 
"O'' = for sequential output (to be written to) 
"I" = for sequential input (to be read from) 
¢see Turbo Basic Ref, pg 282, for other choices| 


1 
q 
t 
{ 
FileNum% 
q 
q 
q 


= file number to be used in read/write calls 
FileNameS = Drive: /path/name. ext of file 
RecLength% = record size in bytes, 1 to 32767, typical = i128 


SUB OpenFile(modeS$ ,FileNum®%,FileNameS ,RecLength*%) 


LOCAL errnum%, ans$ 
t 


Opensequence: 


ON ERROR GoTo ErrorHandling 

OPEN modeS ,FileNum%,FileNameS ,RecLength% 
ON ERROR GoTo 0 

EXIT SUB 


ErrorHandling: 


Beep 2 

errnum% = ERR 

Pelect Gase errnums 
Base 752, /5,/76 


Print "FILE or DRIVE PATH Not Found for ... ", FileNames 


Input "Re-Enter Drive: path name. ext ----- > "; FileNames 
Case 68,71 

Print "DISK DRIVE NOT READY for ", FileName$ 
Case Else 


Print "Opening File Gave Error Number ", errnum% 


End Select 
Input "(R)etry or (H)alt Program "; ans$ 
If UcaseS(ans$) = "H', Then 

Resume HaltSequence 


Else 
Resume OpenSequence 
mad If 
HaltSequence: ‘terminates main program 
t 


Print "Program Halted!!!" 
Print "Subroutine OpenFile with Error Number 
Print '' encountered at program address '', ERADR 
END 

END SUB 


4 LJ * . 
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FO EN GN ER ON GS GN GY EN EN EN ERED EDEL ED ED ED ER EN ER ELEC EV ELECTED ED EC ED COED OU EV AR EV RAD ANED EN EN GN AN EN ON BD 


'LU Back Substitution 
t 


t 


| Ref: Numerical Recipies, by WH Press, page 3/7 
SUB LUBKSB( A#( 2) ,N%, INDX%( 1) , BC 1)) 
LOCAL 11%,1%, j%, sum; 
1i%= 
For it=1 to N% 
11%=Indx%( i%) 
sum=B#( 11%) 
BC 11% )=BECI% 
If ii14<>0, then 
For j%=iir to i-1 
sum#=sumi#-A#Ci%, 1%) BC j%) 
Next j% 
Elself sum#<>0.d0 Then 
Livn=io 
Enid elt 
B#( 1%) =sum# 


{say 
tay 


t 0/ 
seSlLrnums 


aealantwa'nnton'on'antsalsaton’enton's alanlwaton 
FV FR ER FR FN FV EX FR EN GK EXD FU FV FD EDEN E 


t als ataalsn'on'se' on sla nton'e ao ele a one eo le eo no rarer Pro rae rs Ps Ns nt tat tat act at tat le lt ew oe ees wee tet el let la tee wlan as tm Ps es oles it al wes eae wesc en wale a salon wale wl, 
ca AV EN AVL AR EL SV FD CR PH EL FL EN GDR FX IR EL Hd ED FY ray @yx ev er ee rer eri) @yeu ay ee rier ky FV EV IL EV GR EV GL ED EL FEL ER ER ED FE EL GN AL ELEN YN EV ED Ed rd Fd 
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Next i% 
? 


For i%=N% to 1 Step -1 
sum#=B( 1%) : 
If i%<N%, Then 
For j%=ixtl to N% 
sum#=sum#-A#(i%, j%) *Bi#C 4%) 
Next j% 
Jeipuel Ile 
BHC i%)=sum#/A#( 1%, 1%) 
Next i% 
END SUB 


l eedeledistololideclotectooliivedntotctoledednlokvtcicdededetevoinivivinnicckivdocccdcokclteiocedcdcntete: 


Vostsccctetstoncccccowocsccosccccldkkiccdc“s“dTlecc“csTac nso“ ccc: 


'LU Decomposition Subroutine 
’ 


Ref: Numerical Receipes, by WH Press, pages 34-36 
SUB LUDCMP( A#(2) ,N%, Indx%( pap 

LOCAL vvi##() ,Dum#,i%, j%,k%, imax%, aamax# , sum# 

STATIC nmax%, tiny 

nmax%s=10 

tiny#=1. OD-20 


DIM Vive emmax. 


be 
Form 1 ,—) touNs 
aamax#=0 


For j%=1 to N% 
If ABSCA#(i%, j%) )>aamax#, Then aamax#=ABS(A#(i%, j%)) 
Next j% 
If aamax#=0, Then 


Print Matrix has zero mow meets 
Exes Sub 
IZpqvel lice 
VV#Ci%)= 1/aamax} Save COW Seal faerenr. 


Next i% 
For j=l to N% 
eyo lee tila 
FOr ini) te jo- i 
Sum#=A#( 1%, 4% 
If i°;>1 Then 
For ki=1 to 14-1 
Sum-=Sumi} - Ai (1%, k%) *A#C KS, J%) 
Next k% 
A##( 1%, j%)=Sum# 
Spgvels 1a 
Next i% 
Ends it 
continuing j loop, searching for largest pivot 
aamax#=0 
FOL dio jos coun - 
Sum#=A#( 1%, j%o) 


[si 
ON 


If j%>1, Then 
For k%Z=1 to j%-1 
Sum#=Sum# -Aff(i%,k%) *ACK%, 5%) 
Next ik, 
A#( 1%, 4% )=Sum#} 
End: af 
Dum/=vvi#( i%)*ABS(Sum#) ‘figure of merit for pivot 
If Dum#>aamax#, Then 
imax%=1% 
aamax#=Dum}} 
Pratt 
Next i% 


ee) 4c imax. chen 
For k%=1 to N% 
Dum#=A#}( imax%,k%) ‘interchange rows if reqd 
Ad( imax ,K%)=AH( 4%, K%) 
A#( 4%, k%) = Dum} 


Next k% 
D%=-D% ‘change parity of D% 
vvi#( imax%)=vv#( 4%) 
Bad. it 
Indx%( j%)=imax% "divide by pivot 


If (4%<>N%), Then 
Peer Gee] 0, Then Ag( 145 14 )=tiny? 
Dum=1/A#( j%, j%) 
For i%=jrtl to N% 

A#( 1%, 3% )=AHCi%, 3%) *Dum;? 

Next i% 

ng 1 £ 

Next j% 


If A#(N%,N%)=0, Then A#(N%,N%)=tiny# 
END SUB 
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3, 609863.128, 4076611.345, 33.9, 0.0, 0.25, PACKARD 
emu I 289, 4055915. 264, 137.2, 0.0, 0.25, HAYES 
pee 00434.009, 4051260.014, 10.0, 0.0, 0.25, DOPPLER 
6, 596669.490, 4051826.411, 134.5, 0.0, 0.25, HANK 
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APPENDIX B. GEOMETRICAL ACCURACIES OF THE MONTEREY 
BAY MINIFRANGER NETWORK 


Monterey Bay Vicinity 





Figure 13. STATION COMBINATION TREVOR-WATS-HAYS-HANK 
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Figure 1[4. STATION COMBINATION 
TREVOR-PACKARD-DOPPLER-HANK 
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Figure 13. STATION COMBINATION WATS-PACKARD-HAYS-HANK 
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Figure 16. STATION COMBINATION PACKARD-HAYS-DOPPLER-HANK 
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Figure 17. STATION COMBINATION TREVOR-WATS-HANK 
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Monterey Bay Vicinily 





Figure 18. STATION COMBINATION TREVOR-PACKARD-HANK 
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Monterey Bay Vicinity 





Figure 19. STATION COMBINATION WATS-HAYS-HANK 
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Monterey Bay Vicinity 





Figure 20. STATION COMBINATION TREVOR-PACKARD-HAYS 
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